Hierarchical Spatial Organization of Geographical Networks 
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In this work we propose the use of a hirarchical extension of the polygonality index as a means 
to characterize and model geographical networks: each node is associated with the spatial position 
of the nodes, while the edges of the network are defined by progressive connectivity adjacencies. 
Through the analysis of such networks, while relating its topological and geometrical properties, it 
is possible to obtain important indications about the development dynamics of the networks under 
analysis. The potential of the methodology is illustrated with respect to synthetic geographical 
networks. 

PACS numbers: 



Geographical complex networks have been used to 
model biological systems - e.g. bone structure [l[ and 
mammalian cortical areas Q • In these system, the phys- 
ical proximity and communication between neighboring 
elements are critical to guarantee proper development 
and biological function. Another important property of 
these systems is its spatial organization. It is known that 
many biological systems depend on the proper adjacency 
of cells for correct development - e.g. cell communication, 
distribution of cells in retina, kidney structures, among 
many others. 

Despite the large set of tools and measurements used to 
characterize the networks that underly these systems Q , 
relatively little attention has been given to the spatial 
organization of the nodes. In order to address this issue, 
we propose the use of the polygonality index [1] - a robust 
measurement able to quantify the spatial order of system 
of points - over the nodes of a geographical network. 

The polygonality index assigns to each node of the net- 
work a value that indicates the amount of the organiza- 
tion around the node with respect to its neighbors. This 
value is computed based on the angle on formed between 
successive neighbors of a node, as exemplified in Fig. Q] 
and expressed by the following equation: 
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where k is the node under analysis, Nk is the number 
of neighbors of the node k and (3k is a parameter whose 
value can be fixed - to characterize specific spatial ar- 
rangements - or can vary accordingly to the number of 
the neighbors of the nodes, i.e, j3 = 2ir/Nk- A fixed value 
of f3 allows one to identify whether the spatial position of 
the adjacent nodes obeys a specific arrangement (hexag- 
onal, orthogonal, etc.). For example, if j3 = tt/3, it is 
possible to quantify how much the nodes of network dif- 
fer from an hexagonal arrangement. When the value of 
[3 depends on the number of the neighbors, the polyg- 
onality index indicates if the angles between neighbors 
are equally distributed. In both cases, the polygonality 
index varies between (total lack of spatial order) and 1 



FIG. 1: Quantifying the spatial organization around a 
node. All the angles on formed between successive neighbors 
(qi) of the node pt are used to compute the polygonality index 
(see Eq. 1). 



(fully organized system). 

In this article we extend the concept of polygonality 
in order to allow the quantification of the spatial order 
in a complex network including several hierarchical lev- 
els, allowing the quantification of the spatial organization 
along a wider neighborhood around each node. Figure [2] 
illustrates the computation of such a polygonality index 
for three different hierarchical levels of a node. Note that 
virtual edges^ are established between the central node 
and its hierarchical neighbors. The angles between suc- 
cessive virutal edges are taken into account to compute 
the polygonality index, as expressed in Eq. 1. 

Figure [3] shows six networks and respective polygonal- 
ity indices. In (a-c) the networks are perfectly arranged 
in hexagonal lattices, while in (d-f) the networks origi- 
nated from a hexagonal lattice, but the positions of their 
nodes were perturbed in order to reduce the overall order. 
In (a) and (d) the first hierarchal level was considered to 
estimate the polygonality index, while in (b) and (e) the 
second level was used. The results in (c) and (f)consider 
the third level. In all these images the nodes were colored 
accordingly to their polygonality index, as expressed in 
the color-bar. 

Analyzing the values of the polygonality index in Fig- 
ure [31 it is possible to note that in (a), as the nodes and 
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FIG. 2: Hierarchical neighbors used to estimate the spatial organization. Hierarchical neighbors 1 (a), 2 (b) and 3 
(c). The virtual edges (in blue) between the central node and its hierarchical neighbors are used to estimate the angles en for 
polygonality index computation (see Eq. 1) . 



its neighbors are arranged in a hexagonal way, the polyg- 
onality index for all nodes is the maximum possible. In 
(b), the angles defined between successive neighbors of 
the nodes are constant and equal — tt/6. As conse- 
quence, the polygonality index is also at its maximum 
value. On the other hand, in (c) the angles between 
successive neighbors of the third hierarchy are no longer 
constant, but the polygonality index is the same for all 
nodes. In the case of the perturbed lattices (Figure [3]^ d- 
f)), the polygonality index for all hierarchies are, as ex- 
pected, lower than the regular lattices. 

All in all, the presented measurements allowed a com- 
prehensive characterization of the spatial organization of 



geographical complex networks. They can also be com- 
pared with several other measurements related to the dy- 
namics and connectivity of the network, paving the way 
to obtaining new insights about the topological and geo- 
metrical properties of these networks. 
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FIG. 3: Hierarchical polygonality index for a perfect hexagonal lattice (a-c) and for a perturbed hexagonal lattice (d-e). The 
hierarchial level (R), the mean (fi) and standard deviation (a) values of the polygonality index are: (a) R = 1, fi = 1, a = 0; 
(b) J? = 2,/i = 1, a = 0;(c) R = 3, fM = 0.73, a = 0; (d) _R = 1, /x = 0.63, a = 0.07; (e) i? = 2, /i = 0.64, a = 0.05; (f) i? = 3, 
fj, = 0.61, a = 0.05. 



